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We apply the Skyrme model to dense hadronic matter, which provides a unified approach to high density, valid 
in the large Nc limit. In our picture, dense hadronic matter is described by the classical soliton configuration with 
minimum energy for the given baryon number density. By incorporating the meson fiuctuations on such ground 
state we obtain an effective Lagrangian for meson dynamics in a dense medium. Our starting point has been the 
Skyrme model defined in terms of pious, thereafter we have extended and improved the model by incorporating 
other degrees of freedom such as dilaton, kaons and vector mesons. 



1. Introduction 

At high temperature and/or density, hadrons 
are expected to possess properties that are very 
different from those at normal conditions. Un- 
derstanding the properties of hadrons in such ex- 
treme conditions is currently an important issue 
not only in nuclear and particle physics but also 
in many other related fields such as astrophysics. 
Data from the high energy heavy ion colliders, 
astronomical observations on compact stars and 
some theoretical developments have shown that 
the phase diagram of hadronic matter is far richer 
and more interesting than initially expected. Lat- 
tice QCD calculations have been carried out suc- 
cessfully at high temperature, however similar 
calculations at high density have not yet been 
possible. Theoretical developments have unveiled 
such interesting QCD phases as color supercon- 
ductivity. Moreover effective theories can be de- 
rived for these extreme conditions, using macro- 
scopic degrees of freedom, by matching them to 
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QCD at a scale close to the chiral scale A^ ^ 
AwU - 1 GeV. 

We have followed a different path to dense 
matter studies 1 2'3'4'5'6't by using as our start- 
ing point a model Lagrangian, in the spirit of 
Skyrme, which describes hadronic matter and me- 
son dynamics respecting the symmetries of QCD. 
The parameters of the model are fixed by me- 
son dynamics at zero baryon number density. A 
la Skyrme, baryons arise from a soliton solu- 
tion, the skyrmion, with the topological wind- 
ing number describing the baryon number. In 
our scheme dense matter is approximated by a 
system of skyrmions with a given baryon num- 
ber density whose ground state arises as a crys- 
tal configuration TT. Starting from this ground 
state our approach provides insight on the intrin- 
sic in-medium dependence of meson dynamics. 
We have studied (i) the in-medium properties of 
the mesons and (ii) the role of the other degrees of 
freedom besides pions in the description of matter 
as it becomes denser. 
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2. Model Lagrangian 

The original Skyrme model Lagrangian[7] reads 



f2 1 

4 J - 32e^ - (1) 



where = C/^(9pt/ and [/ = exp(if • tt) G S'[/(2) 
is a nonlinear realization of the pion fields, f^^ 
the decay constant and TOtt the pion mass. The 
second term with e, the Skyrme parameter, was 
introduced to stabilize the soliton solution. 

The dilaton field x can be incorporated in 
the model to make it consistent with the scale 
anomaly of QCD 8 9 . The Lagrangian then 
becomes modified as 
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Note the different powers of {x/ fx) front 
of each term. The last line is the La- 
grangian for the free dilaton field, where V{x) — 

{mlflimx/k)\Hx/fx) - i) - i), rn^ is the 
dilaton mass and its decay constant. 

The vector mesons, p and w, can be included 
into the model as dynamical gauge bosons of a 
hidden local gauge symmetry which requires the 
doubling of the degrees of freedom as U = ClCfl- 
One of such Lagrangian is, for example JU' , 
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di,Pfj, + gPi_i X pi,, iVf^ty = d^uji, - d„uj^, and where 
is the topological baryon number current. The 



quartic Skyrme term of ^ is not present, because 
its stabilizing role is played here by the vector 
mesons. 

One may decouple the vector mesons and/or 
dilaton field from the other by making the corre- 
sponding particles infinitely heavy. In this limit, 
the dilaton field is frozen to its vacuum value 
X — fx ^^'^ '^ho mesons are constrained to 
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as kind of composite particles made of two pions. 
Then, the model Lagrangian (0) reduces to (O 
or where the kinetic energy term of the rho 
vectors becomes the Skyrme term. 

3. Dense Skyrmion Matter 

These nonlinear meson Lagrangian supports 
soliton solutions, skyrmions, carrying nontrivial 
topological winding numbers. Once we accept 
Skyrme's conjecture of interpreting the winding 
number as the baryon number, we may describe a 
dense baryonic matter as a system made of many 
skyrmions. 

In the classical picture, the lowest energy state 
of the multi-skyrmion system is a crystal and 
there has been intensive work in late 80's with the 
model containing pions only m. Well-separated 
two skyrmions have lowest energy when they are 
relatively rotated in the isospin space about an 
axis perpendicular to the line joining their cen- 
ters. We can imagine that the lowest energy 
state of skyrmion matter in a relatively low den- 
sity is in an FCC (face centered cubic) crystal 
phase where well localized single skyrmions are 
arranged on each lattice cite in a way that 12 
nearest skyrmions have such lowest energy rela- 
tive orientations. At higher density, on the other 
hand, skyrmion tails start overlapping and the 
system undergoes a phase transition to a more 
symmetric configuration the so called the "half- 
skyrmion" cubic crystal. There, one half of the 
baryon number carried by the single skyrmion is 
concentrated at original FCC site where U = —1 
while the other is concentrated on the links where 
U = +1. Now, the system has an additional sym- 
metry with respect to U ^ — ?7, which results in 
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Figure 1. Average values of a and x/ fx the 
lowest energy crystal configuration at a given 
baryon number density. 



the vanishing of {(j){U = a + if ■ ip). In the lit- 
erature, it is often interpreted as the restoration 
of the chiral symmetry. However, it is only the 
average value of a over the space that vanishes, 
while the chiral circle still has a fixed radius /tt. 
We call this phase as "pseudogap" phase to dis- 
tinguish it from the genuine chiral symmetry re- 
stored phase, where the chiral circle shrinks to a 
point at cr = 0, = 0. 

In the model with dilaton field which plays 
more or less a role of the "radial" field for U, 
the restriction on chiral radius becomes released. 
The pseudogap phase still remains as a tran- 
sient process unless the dilaton mass is sufficiently 
small^. Shown in Fig.l is a typical numerical 
results on the average values of a and x/ fx- 
the figure, the explicit scale of the baryon num- 
ber density, which depends strongly on the model 
parameters, is not shown. One can see that, as 
the density increases, the average value of a drops 
quickly and reaches zero at the density pp, where 
the system changes to the pseudogap phase. As 
we increase the density further, the system re- 
mains in the pseudogap phase for a while but the 
average value of x/ fx keeps going down slowly. 
At the density pc, the {xl fx) ^ phase and the 
{xl fx) — phase have the same energy. Then, at 
higher density than pc, the latter comes to have 
lower energy and finally chiral symmetry is re- 
stored. 



The vector mesons (especially the omega me- 
son) also play important roles in dense skyrmion 
matter. Numerical results obtained in various 
models are presented in Figure 2. In the 7rp 
model, as the density of the system increases (and 
the lattice constant decreases) , the system energy 
per baryon, E/B, changes slightly. Its value is 
close to the energy of a single skyrmion up to 
somewhat high densities. The main reason for 
this is that the size of the skyrmion is very small 
because of the absence of strong repulsive terms 
in the model. Thus, the skyrmions in the lattice 
interact only at very high densities where their 
tails overlap. In the i^px model without oj, the 
dilaton field plays an important role. Skyrmion 
matter undergoes an abrupt phase transition at 
the density where the expectation value of the 
dilaton field vanishes (x) = 0. 

In the npujx model, the situation changes dra- 
matically. Above all, to provides a strong repul- 
sion which inflates each single skyrmion. The tails 
of the bigger skyrmions overlap providing attrac- 
tion to the system in the intermediate range. 

In both the Trpcj and the Trpujx models, at high 
density, the interaction reduces E/B to 85% of 
the B = 1 skyrmion mass. This value should 
be compared with 94% in the irp model. In the 
Trpx- model, E/B goes down to 74% of the B = 1 
skyrmion mass, but in this case it is due to the 
dramatic behavior of the dilaton field. On the 
other hand, the omega suppresses the role of the 
dilaton field. It could provide only a small at- 
traction at intermediate densities. Moreover, the 
phase transition towards its vanishing expecta- 
tion value, (x) — does not take place. Instead, 
its value grows at high density! 

The reason for this can be found in the role 
played by omega in (jSJ. In the static configura- 
tion, omega produces a potential, whose source is 
the baryon number density, which mediates the 
self-interaction energy of the baryon number dis- 
tribution. Thus, unless it is screened properly by 
the omega mass, the periodic source filling infi- 
nite space will lead to an infinite self-energy. To 
reduce the energy of the system, the effective lu 
mass must grow at high density, for which x must 
grow too. Note the factor (x//x)^ the omega 
mass term in Lagrangian Q. 



<x'f,> 1/ 

\ 




j\ 





E/B 


\. 






\ jtpD»- model ^- 








0.5 




npoiz- model 





1.0 2.0 0.0 1.0 

FCC latlice constant FCC lattice constant 



Figure 2. Numerical results in the model with 
vector mesons and dilaton field. 



On the other hand, these classical crystalline 
structures are quite far from normal nuclear mat- 
ter which is known to be a Fermi liquid at low 
temperature. In order for skyrmion matter to be 
identified with nuclear matter we have to quan- 
tize and thermalize the classical system. Since it 
is a system of solitons in a meson field theory, 
it is not sufficient to quantize the meson fiuctua- 
tions. We need to introduce and quantize proper 
collective variables not only to complement the 
broken symmetries of the whole skyrmion system 
but also to describe the dynamics of the single 
skyrmions in order to obtain a realistic picture 
of nuclear matter. For extended objects, we may 
need an infinitely large number of dynamical vari- 
ables such as the positions of their center of mass, 
their relative orientations, their sizes, their de- 
formations, etc. Among them, those degrees of 
freedom that describe translations and rotations 
of the single skyrmion play the dominant role at 
low energy. Thus, we need at least 6 variables 
for each skyrmion. For a multi-skyrmion system, 
the simplest way of introducing collective vari- 
ables for the position and orientation of each sin- 
gle skyrmion is through the use of the product 
Ansatz, the old idea of SkyrmejZj- There, a multi- 
skyrmion solution can be approximated by prod- 
ucts of single skyrmion solutions centered at the 
corresponding positions and rotated to have the 



corresponding orientation. However, the product 
Ansatz works well only when the skyrmions are 
sufficiently separated. Furthermore, due to the 
non-commutativity of the matrix products, it is 
difficult to use it in multi-skyrmion systems. 

Another scheme which can be used to study 
multi-skyrmion systems is the Antiyah-Manton 
Ansatz |12| . In this scheme, skyrmions of 
baryon number N are obtained by calculating 
the holonomy of Yang-Mills instantons of charge 
N ^ which has been used in describing success- 
fully few-nuclcon systems and also nuclear mat- 
ter. One advantage of the Atiyah-Manton ansatz 
is that it provides a natural framework to in- 
troduce the proper dynamical variables for the 
skyrmions through the parameters describing the 
multi-instanton configuration. Contrary to the 
non-commutative product ansatz, some multi- 
instanton solutions are given in a commutative 
manner. Furthermore, multi-instanton solutions 
have been investigated widely and many useful 
solutions have been found. 

We have tried this idea in Ref.J^. 'tHooft 
ansatz on multi-instanton solution is adopted to 
produce the FCC instanton crystal. To avoid the 
divergence coming from the infinite number of 
instantons and to incorporate the relative orien- 
tations to each instantons, the ansatz is slightly 
modified. Then, the Antiyah-Manton procedure 
is applied to transform the instanton crystal into 
a skyrme crystal. As expected, the resulting 
skyrmion crystal is FCC for low baryon number 
density. At high density, however, it becomes a 
half- skyrmion CC (approximately). 

One great advantage of this procedure to gen- 
erate the skyrmion crystal is that it is really made 
up of single objects located at specified positions 
and with specified rotations. Now, we can, for 
example, move a single skyrmion and investigate 
how the system changes. Shown in Fig. 3 is the 
potential energy V{d) when a single skyrmion is 
moved away by a distance d from its stable posi- 
tion. Two extreme cases are shown. In the case of 
a dense system {Lp(j(j = 5.0), the energy changes 
abruptly. For small d, it is almost quadratic in d. 
It implies that the dense system is in the crys- 
tal phase. On the other hand, in the case of di- 
lute system {LpQQ — 10.0), the system energy 
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Figure 3. The energy cost V{d) to shift a single 
skyrmion ft'om its stable position by an amount 
d along the edge of FCC lattice. 



remains almost constant up to some distance d, 
which implies that the system is in a gas or liq- 
uid phase. If we allow all the variables to vary 
freely to seek the minimum energy configuration, 
the system will end up in a disordered phase, in 
which a few skyrmions will form a finite cluster. 
Furthermore, we will be able to investigate the 
thermal properties of the skyrmion system and 
quantize those collective variables. 

4. Mesons in Dense Matter 

By using the same model Lagrangian defined at 
zero baryon number density, we may investigate 
the properties of the mesons in dense matter. The 
Lagrangians CQ)-® are originally constructed to 
describe the meson dynamics in vacuum, which is 
defined as 

U=l, X = fx and p^=cj^ = 0. (5) 

The fluctuations on top of this vacuum can de- 
scribe the corresponding particles in zero baryon 
number space, for which we determine the physi- 
cal parameters of the model Lagrangian, such as 
the particle masses, decay constants and coupling 
constants. 

Now, as we have described in the previous sec- 
tion, we have another classical background config- 
uration for the matter with non vanishing baryon 
number density. Let's denote the solution as U,^-., 



Xj-g-j , P^'q-, and w^j,^ , respectively. We can incorpo- 
rate the meson fluctuation on top of this back- 
ground as 

U = yC^C/(0) V^, X = X(0) + X* , 

where the 'starred' flelds are describing the cor- 
responding particles in medium. 

We will illustrate the basic strategy of the ap- 
proach by using the simplest model with only 
pious. Substituting eqs.® into the Lagrangian 
density ^ and keeping the terms up to the sec- 
ond order in the fluctuating flelds, we obtain 

where Gabi^") and cr(r) are the background poten- 
tials provided by the dense skyrmion matter. For 
the fluctuating pious in free space, both poten- 
tials are just 1. This Lagrangian tells us how the 
properties of the pious change in the dense mat- 
ter. For example, Gabi?) in front of the pion ki- 
netic term can be absorbed into the wavefunction 
renormalization, which can be interpreted as a ra- 
tio between the local effective pion decay constant 
in dense medium to that in free space. The pion 
mass term gets a similar correction from dense 
matter. To the first order in the background po- 
tentials, we can translate iQ into an effective La- 
grangian for the pious in dense medium with in- 
medium physical parameters as 

In Figure 4 we show the ratios of the in-medium 
parameters of pious relative to its free-space val- 
ues. The pion mass is almost constant at low 
density but decreases as the density increases to 
much higher values. The pion decay constant 
drops quite fast but after some density it stays 
with value '--^ 2/3 indicating that the system is in 
the pseudogap phase mentioned in Section 3. A 
similar process can be applied after incorporat- 
ing the dilatonjSj, where the vanishing of x can 
restore the chiral symmetry completely. 

One may treat the interaction of the fluctuating 
flelds with the background potentials in a more 
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Figure 4. In medium properties of the pions as a 
function of the baryon number density. 



systematic way. As an example, in Ref. X the in- 
medium pion velocity is studied, where the back- 
ground interactions are taken into account up to 
the second order. The breakdown of the Lorentz 
symmetry due to the presence of medium makes 
the pion velocity deviate from that in free space. 

5. Summary 

We have developed a unified approach to dense 
matter within the Skyrme philosophy, where sys- 
tems of baryons and mesons can be described by 
a single Lagrangian. In our approach dense bary- 
onic matter is approximated by skyrmion mat- 
ter in the lowest energy configuration for a given 
baryon number density. By incorporating in it 
fluctuating mesons we can get some insight on 
meson dynamics in a dense medium. Our ap- 
proach enables us to study this dynamics beyond 
the first order in the baryon number density. One 
can continue to work in this direction by incor- 
porating more degrees of freedom, by improving 
the way of treating matter beyond the crystal so- 
lution, and so on. 

However, before closing the presentation, we 
must clearly lay down the scope of our work. We 
do not claim that the results obtained at present 
describe reality. The most fundamental problem 
we phase is that our "ground state" for matter is 
a crystal not a Fermi liquid. Our aim has been 
to assume a state for matter, given by a classi- 



cal solution of a theory considered to be valid at 
large Nc, and have studied the implications for 
its excitations. Our work should be taken as rep- 
resenting the first step towards a more realistic 
treatment of a dense matter theory. 
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